UG-FT- 75/97 
hep-ph/9708395 



m u + m d 



from QCD-Hadron Duality!! 



Joaquim Prades 

Departamento de Ffsica Teorica y del Cosmos, Universidad de Granada 
Campus de Fuente Nueva, E-18002 Granada, Spain. 



Abstract 

The value of the light quark masses combination m u + is analized 
using QCD-Hadron Duality. A detailed analysis of both the perturbative 
QCD [to four-loops] and the hadronic parametrization needed is done. The 
result we get is [m u +m d )(l GeV 2 ) = (12.8±2.5) MeV ([m u +m d ]{A GeV 2 ) = 
(9.8 ± 1.9) MeV ) in the MS scheme. 
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We are interested in the value of the light quark masses that appear in the 
Standard Model. These are the QCD Lagrangian quark masses or Current Al- 
gebra quark masses. This work is an update of |1J putting some emphasis in the 
QCD-Hadron Duality and including the recently calculated a 3 QCD corrections 
[T2J. A collection of the relevant references to earlier work can be found in [I]. 

The QCD Sum Rule (SR) technique has been discussed several times in 
this conference. (See the talks by Narison and Jamin for instance.) The rele- 
vant two-point function for calculating m u + rrid is the correlator of two diver- 



gences of the axia 

i G (j p ) = i-(o-) 



•vector current with the quantum numbers of the pion field 
(See also the talk by de Rafael.) 
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x) = [m u + nid] [ui'fsd]. The QCD Sum Rule tech- 
nique exploits the analytic structure of ty 5 (q 2 ) combined with Cauchy's theorem. 
The specific type of SR we shall use is the so-called Finite Energy Sum Rule 
(FESR). For the discussion of quark masses using QCD-Hadron Duality there 
are two relevant FESRs. These are the so-called Oth-moment 
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and lst-moment 
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respectively. The upper limit of the integral, s, is the onset of the QCD contin- 
uum. The hypothesis of QCD-Hadron Duality that we have used to obtain (fj) 
and (H) can be expressed as follows. One expects that there exists some interme- 
diate region [s belongs to that region], such that above this region the description 
of inclusive enough quantities like ^^(q 2 ) are well approximated by perturbative 
QCD, i.e. local duality is approximately satisfied. In the intermediate region one 
can just expect to have duality with perturbative QCD for suitable averages of 
inclusive quantities like the moments in (H) and @, i.e. global duality. Below 
that region one doesn't expect any kind of duality. 

Let me discuss the inputs needed in these FESRs. The rhs in @ and (H) are 
given by 
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with n — 1 and n = 2 respectively, and \t) = 4j5^s{t)- For large enough 
Q 2 = — g 2 = — £, ^^(Q 2 ) is known in QCD up to four-loops [order a^j from 
[0, 0]. Once the global light quark masses factor scaled at s together with the 
energy dependence and normalization factors have been pulled out, what remain 
from (H) is the QCD perturbative series 1 + Ri(s). This is supplemented by power 
corrections a la SVZ ||] parametrized by Wilson coefficients Ci and condensates 
(Oi). The values of the power corrections coefficients are C 4 (0 4 ) = (0.08 ± 0.04) 
GeV 4 and C 6 (O e ) = (0.04 ± 0.03) GeV 6 Q. 

To study the convergence of the truncated perturbative series, we have used 
three different ways of resumming the QCD series. The first one, that we shall 
call Perturbative, is obtained resumming the renormalization group logs after the 
Cauchy's integral in (|]) is done. In the second one, which we shall call Improved, 
the large running of the coupling constant around the Cauchy's circuit (from 
to 2tt) is taken into account by resumming the renormalization group logs 
before the Cauchy's integral in (|]) is done ||. The third way of resumming is 
obtained by applying the Principle of Minimal Sensitivity (PMS) |§ to the QCD 
expression for ^/^(Q 2 ) before the Cauchy's integral is done. We have applied 
the renormalization group running to a s (m T ) = 0.35 ± 0.02 at three-loops to 
obtain a s (s) around the Cauchy's circuit. Below we show the QCD perturbative 
series 1 + R2{s) for s = 2 GeV 2 obtained using the three ways of resumming. 



Perturbative 
Improved 
PMS 



1 + 0.60 + 0.38 + 0.19 + •• • 
1 + 0.68 + 0.20 + 0.07 + ••• 
1 + 0.85 + 0.12- 0.01 + ■•• . (5) 



The second term in each line is the order a s correction, the third is the order a 2 
and the fourth the order correction. The comparison of the series in @ shows 
that PMS is the best behaved but also shows how the three series show good 
convergence if one considers corrections up to order a^. The numerical difference 
between them is less than 10 % at that order. 

Let me now discuss the inputs needed for the lhs of (Q) and (|[). The integrand 
lm^5(t) is an inclusive cross-section which, ideally, should be determined in terms 
of the hadronic states that couple to the vacuum through the divergence of the 
axial- vector current d^ 1 A l ~ l2 {x), i.e. the 7r-pole, 37r-continuum, 57r-continuum, 
7r'-resonances, NN , ■ ■ ■. Apart of the pion pole, we have not complete access to 
this spectral function yet. We do however have several pieces of information on 
it. In the low energy end, Chiral Perturbation Theory (CHPT) predicts 

-Im* 5 (i) = 2f 2 m*5(t-ml) 
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where f w = 92.4 MeV in this normalization and the dots stand for higher thresh- 
olds contributions. The first term is due to the pion pole and is completely fixed 
by the pion mass and its decay coupling constant /„■, this will become important 
afterwards. When t — > 0, the spectral function p^(t) can be obtained at low- 
est order in CHPT. At lowest order in Generalized CHPT, one gets the CHPT 
result times some factor, see Stern's talk. In fact, no symmetry allows to pre- 
dict p^(t = 0) and it receives corrections to all orders in CHPT proportional 
to pion and kaon masses. We shall let this normalization factor free and require 
QCD-Hadron Duality to fix it, see below. 

For values of t above a few hundred MeVs, CHPT is not anymore at work and 
we enter the resonance region. Here we have some experimental information to 
help us, for instance we know that there are resonance states with the quantum 
numbers of the pion which couple to the 37r-continuum, 57r-continuum, • • •. There 
are also two-body thresholds like pir, p'n, pir', NN, and so on. The contributions 
of the resonances are in fact leading in the 1/N C expansion [p], |J and have to be 
included if we want to obtain QCD-Hadron Duality. In the Review of Particle 
Physics 0, we can find two resonances with the pion quantum numbers. Their 
masses and widths are: M 1 = [1300 ± 100] MeV, T t = [400 ± 200] MeV, M 2 = 
[1770 ± 30] MeV, T 2 = [310 ± 50] MeV. The most recent experimental data 
where both tt'- resonances are seen is in ||10|| . One can find there the energy 
dependence of the I G (J P ) = 1~(0~) amplitude in the three-pion channel. This 
amplitude is proportional [in the three-pion channel] to the spectral function 
with the quantum numbers of the divergence of the axial-vector current we are 
studying. Unfortunately, the normalization of the coupling of this current to iViV 
is unknown. That data allowed to obtain the masses and widths of the two n'- 
resonances since the non-resonant background in the data was found to be very 
small [ lOfl . The values obtained there agree quite well with the ones reported in 



the Review of Particle Physics ||. Therefore, we shall take the I G (J P ) = l - (0~) 
data in [[Kj as a good estimate of the resonant shape allowing a free relative 
strength between both tt' resonances. 

Another piece of information to add to our study of the spectral function 
is the success of the resonance dominance in the vector channel (VMD). One 
can obtain the successful VMD predictions using a spectral function with Breit- 
Wigner shapes for the resonances normalized to lowest order CHPT prediction [p], 
||. The difference in this case with respect to the divergence of the axial- vector 
channel one is that the normalization at t = of the vector form factor is fixed 
since the vector current couples to the electromagnetic current. The proposal 
we make is to use scalar meson dominance for the I (J p ) = 1~(0 - ) channel, so 
that the spectral function for the three-pion channel is the one in (|6|) with the 
substitution 

p^it)-, p^{t)=Ap^{t) \T{M 1 ,T 1 ,M 2 ,T 2 ;t)\ 2 (7) 
where A is a free normalization factor, p^it) is the lowest order CHPT result 
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and 



J^(M 1 ,T 1 ,M 2 ,r 2 ;t) 
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-Mf + ifiMi -Mg + iT 2 M 2 



Here £ is a complex parameter which fixes the relative strength between the two 
7r' resonances and a possible different absorption by the nuclear target. We also 
included the corresponding Breit-Wigner resonant shape for the p-meson in the 
intermediate pn two-body subchannel which contribution is leading in 1/N C and 
phase space [p], ||. The £ parameter can be obtained from a fit of the I G (J P ) = 
1~(0~) data in |H|. We get a good fit to the shape of the I G (J P ) = 1~(0~) data 
by fixing f = 0.234 + i 0.1. 

Let me now apply the technique to obtain m u + using all the information 
above. We only use the three-pion channel. The onset of the QCD continuum s 
in (0) and (|3|) is fixed by demanding a good QCD-Hadron Duality between 
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The plots of the ratios in are shown in Figure |I|. The curves labeled as 
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Figure 1: Duality Ratios in (|9|). See text for explanation of the curves. 



Perturbative, Improved, and PMS correspond to 7£qcd(s) obtained using the 



different ways of resumming the QCD series explained above. Observe the small 
difference between the three curves. The effect of varying a s and the condensates 
in the range quoted above is less than 5 %. Other non-perturbative contributions 
like for instance the possible large contribution of instantons mostly cancel in 
IZqcb(s). In particular, the result in [12|] produces less than 5 % change in 



7£qcd(s). After the good control we have on the QCD counterpart, a violation 
of global QCD-Hadron Duality larger than say 10 % is unexpected and certainly 
would have to be understood. 

Notice that duality cannot diiferentiate between shapes of spectral functions 
which have the area below. On the contrary, the interference of the pion pole 
contribution in @ and the three-pion contribution makes the ratio 7£had(s) sen- 
sitive to changing the global normalization factor A. The normalization factor 
A should be around 1.5 if QCD-Hadron Duality is wanted (We have used the 
PMS resummed series as the best behaved.). In Figure |IJ we show three different 
normalizations: namely, the lower curve is for A = 1 in analogy with the normal- 
ization of the conserved vector current; the curve with normalization factor A = 
1.5 is the one that has the best duality with the PMS resummed QCD series, 
and the higher curve is for A — 5 where more than 20 % violation of duality 
is obtained. Notice also that one cannot expect duality for s larger than [3 ~ 
3.5] GeV 2 since from these energies on there are two-body thresholds like NN, 
pn' [which couples dominantly to 57r-continuum], • • •. None of these intermediate 
states have been included in our parametrization or in the one in ||11|| . These 



contributions are leading in the 1/N C counting and dominant by phase space. 



This was overlooked in |TTJ where duality in the three-7r channel was required at 
energies much higher than 4 GeV 2 . More on this issue in ||. 

Once we have fixed the duality region from Figure |1] to be around 2 GeV 2 we 
can resolve for the light quark masses in that region. Using the 0-th moment (Q), 
we get the masses in Figure 0. There we plot the MS m u + m d masses scaled to 
1 GeV obtained using the hadronic parametrization with A — 1.5 and the PMS 
resummed QCD series. The value of m u + in the duality region can be read 
off from Figure 0. The result we get is 

[m u + m d ](lGeV 2 ) = (12.8 ± 2.5) MeV (10) 

i.e. 

[m u + m d ](4GeV 2 ) = (9.8 ± 1.9) MeV (11) 

in the MS scheme at order a^. This result updates the one in Jl|. The central 
value is for A = 1.5 and using the PMS resummed QCD series. The error includes 
the hadronic uncertainty allowing for a violation of duality of more than 10 % 
at s = 2 GeV 2 and by varying a s and the coefficients of the power corrections in 
the ranges given. All these uncertainties have been added in quadrature. Notice 
also that the results one gets using the Improved QCD resummed series or the 
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Figure 2: m u + MS masses scaled at 1 GeV. 



PMS one only differ by 0.1 MeV. This shows the good control we have on the 
QCD series. Even though our spectral function is obviously not complete [see 
comments above], one can see from Figure ^| that the value of the masses is quite 
stable between s = 2 GeV 2 and s = 4 GeV 2 . Instanton contributions to the 0-th 
moment for s around 3 GeV 2 are negligible |12[ and have been included in the 
quoted error. Laplace Sum Rules give compatible results. The results in (|10|) 
and ( |TTD are in perfect agreement with the bounds obtained recently in [13 . 



Given the good behaviour shown by the QCD counterpart, further improve- 
ments reducing the error bars in the determination of m u + have to come 
from the very difficult measurement of the I G (J P ) = 1~(0~) spectral function 
normalization. As we tried to show, our present understanding already gives tight 
bounds on the values of m u + m^. 

I would like to thank Hans Bijnens and Eduardo de Rafael for reading the 
manuscript and useful comments. 
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